First-principles studies of ground- and excited-state properties of MgO, ZnO, and CdO polymorphs 
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An ab initio pseudopotential method based on density functional theory, generalized gradient corrections to 
exchange and correlation, and projector-augmented waves is used to investigate structural, energetical, elec- 
tronic and optical properties of MgO, ZnO, and CdO in rocksalt, cesium chloride, zinc blende, and wurtzite 
structure. In the case of MgO we also examine the nickel arsenide structure and a graphitic phase. The stability 
of the ground-state phases rocksalt (MgO, CdO) and wurtzite (ZnO) against hydrostatic pressure and biaxial 
strain is studied. We also present the band structures of all polymorphs as well as the accompanying dielectric 
functions. We discuss the physical reasons for the anomalous chemical trend of the ground-state geometry and 
the fundamental gap with the size of the group-II cation in the oxide. The role of the shallow Zn3rf and Cd4d 
electrons is critically examined. 

PACS numbers: 61.66.Fn, 71.20.Nr, 78.40.Fy 



I. INTRODUCTION 

Currently the optical properties of wide band-gap semicon- 
ductors such as zinc oxide (ZnO) are of tremendously increas- 
ing interest, in response to the industrial demand for optoelec- 
tronic devices operating in the deep blue or ultraviolet spec- 
tral region. Under ambient conditions ZnO is crystallizing in 
wurtzite (w) structure. Its tendency to be grown as fairly high 
residual n-type material illustrates the difficulty to achieve its 
p-type doping. Nevertheless, the potential of ZnO for opto- 
electronics yj but also for spintronics |_2J (e.g. in combina- 
tion with MnO) renders it among the most fascinating semi- 
conductors now and of the near future. 

Heterostructures with other materials are most important 
for optoelectronic applications jj, |j]. One important class 
of crystals for heterostructures with ZnO could be the other 
group-II oxides and alloys with these compounds. Another 
IIB oxide is CdO which however has a much smaller funda- 
mental energy gap 1 5 ] . On the other hand, the group-IIA oxide 
with the smallest cation, MgO, possesses a much larger energy 
gap 1 6]. Consequently, combinations of these oxides with 
ZnO should lead to type-I heterostructures with ZnO (CdO) 
as the well material and MgO (ZnO) as the barrier material. 
However, there are at least two problems for the preparation of 
heterostructures: (i) The group-II oxides MgO, ZnO, and CdO 
in their ground-states do not represent an isostructural series 
of compounds with a common anion, that means the ground- 
states are given by different crystal structures: the cubic rock- 
salt (MgO, CdO) or the hexagonal wurtzite (ZnO) one |5]. (ii) 
The in-plane lattice constants of two oxides grown along the 
cubic [111] or hexagonal [0001] direction are not matched. 

Among several other properties of the group-II oxides that 
are not well understood, is their behavior under hydrostatic 
and biaxial strain as well as their stability: Which are possi- 
ble high-pressure phases? Are there any strain-induced phase 
transitions? Which role do the c/-electrons play for ZnO 
and CdO? Other open questions concern the influence of the 
atomic geometry on the band structure and, hence, on the ac- 
companying optical properties. Furthermore, the band-gap 
anomaly of CdO has to be clarified: While the band gaps 
of the common-cation systems CdTe, CdSe, CdS, and CdO 



show an increase of the fundamental gap along the decreasing 
anion size until CdS, the value for CdO is smaller than that 
of CdS 1 7]. A similar anomaly occurs along the row ZnTe, 
ZnSe, ZnS, and ZnO fT], however, the gap variation is much 
smaller. Such gap anomalies are also observed for III-V semi- 
conductors with a common cation: one example is the row 
InSb, InAs, InP, and InN |8], where the anomaly for InN has 
been traced back to the extreme energetical lowering of the 
N2.S' orbital with respect to the Sb5.s, As4j', and P3.s levels and 
the small band-gap deformation potential of InN tSI- 

In this paper, we report well-converged ab initio calcula- 
tions of the ground- and excited-state properties of the most 
important polymorphs of the group-II oxides MgO, ZnO, and 
CdO. In Sec. II the computational methods are described. 
Atomic geometries and the energetics are presented in Sec. 
Ill for cubic and hexagonal crystal structures. In Sec. IV we 
discuss the corresponding band structures and the electronic 
dielectric functions. Finally, a brief summary and conclusions 
are given in Sec. V. 



II. COMPUTATIONAL METHODS 

Our calculations are based on the density functional theory 
(DFT) 1 9] in local density approximation (LDA) with gener- 
alized gradient corrections (GGA) llOll according to Perdew 
and Wang (PW91) \IV\. The electron-ion interaction is de- 
scribed by pseudopotentials generated within the projector- 
augmented wave (PAW) scheme I12l[l3lll4ll as implemented 
in the Vienna Ab initio Simulation Package (VASP) 1 15|. The 
PAW method allows for the accurate treatment of the first-row 
element oxygen as well as the Zn3t/ and CdAd electrons at rel- 
atively small plane wave cutoffs. 

For the expansion of the electronic wave functions we use 
plane waves up to kinetic energies of 29.4 Ry (MgO in 
wurtzite, rocksalt, CsCl, NiAs, and graphitic-like structure, 
ZnO in zinc blende and wurtzite structure, and CdO in rock- 
salt and wurtzite structure) and 33.1 Ry (MgO in zinc -blende 
structure, ZnO in CsCl and rocksalt structure, and CdO in 
CsCl and zinc -blende structure), respectively. To obtain con- 
verged results for the external pressures (trace of the stress 



tensor) we increased the plane wave energy-cutoff to 51.5 Ry 
uniformly for all materials. The Brillouin-zone (BZ) inte- 
grations in the electron density and the total energy are re- 
placed by summations over special points of the Monkhorst- 
Pack type 1 161 . We use 8x8x8 meshes for cubic systems and 
12x12x7 for hexagonal polymorphs. 

A first approach to qualitatively reliable band structures 
is using the eigenvalues of the Kohn-Sham (KS) equation 
llOll . They also allow the computation of the electronic 
density of states (DOS). We apply the tetrahedron method 
1 17] to perform the corresponding BZ integration with k- 
space meshes 20x20x20 for cubic crystals or 30x30x18 
for hexagonal structures. For the cubic polymorphs the 
frequency-dependent complex dielectric function £(©) is a 
scalar, but it possesses two independent tensor components 
Exxito) = eyy(ft') and ezz{co) in the cases of hexagonal sys- 
tems. In independent-particle approximation it can be calcu- 
lated from the Ehrenreich-Cohen formula |18]. For the BZ 
integration in this formula we use refined k-point meshes of 
50x50x31 for hexagonal structures and 40 x 40 x 40 for cubic 
crystals. In particular, the frequency region below the first ab- 
sorption peak in the imaginary part depends sensitively on the 
number and distribution of the k points. The resulting spectra 
have been lifetime-broadened by 0.15 eV but are converged 
with respect to the used k-point meshes. 



III. GROUND-STATE PROPERTIES 
A. Equilibrium phases 

For the three oxides under consideration, MgO, ZnO, and 
CdO we study three cubic polymorphs: the Bl rocksalt (rs or 
NaCl) structure with space group Fm3m (O^), the B3 zinc- 



blende (zb or ZnS) structure with space group F43m (T^j 
and the B2 CsCl structure with Pm3m (O/'). In the case of 
the hexagonal crystal system we focus our attention to the B4 
wurtzite (w) structure with space group P63mc (Cg^,) but we 
also investigate the B8i NiAs structure with P63/mmc (Dg/,) 
symmetry and a graphitic-like structure with the same space 
group 1 19] for MgO, which we call h-MgO, according to ] 20]. 
The fourfold-coordinated w and zb structures are polytypes 
with the same local tetrahedral bonding geometry, but they 
differ with respect to the arrangement of the bonding tetra- 
hedrons in [0001] or [111] direction. Their high-pressure 
phases could be the sixfold-coordinated NaCl or eightfold- 
coordinated CsCl structures. In the cubic zb and rs phases the 
cation and anion sublattices are displaced against each other 
by different distances parallel to a body diagonal, (1,1,1 )flo/4 
for zb and (l,l,l)flo/2 for rs, respectively. There are also 
several similarities for the three hexagonal phases wurtzite, 
NiAs and h-MgO. In the NiAs structure the sites of two ions 
are not equivalent. For the ideal ratio c/a = \/8/3 of the 
lattice constants the anions (As) establish a hexagonal close- 
packed (hep) structure, whereas the cations (Ni) form a simple 
hexagonal (sh) structure. Each cation has four nearest anion 
neighbors, whereas each anion has six nearest neighbors, four 



cations and two anions. The latter ones form linear chains par- 
allel to the c-axis. In the case of /z-MgO the cations and anions 
form a graphitic-like structure (2l|j] with Bk BN symmetry. In 
comparison to wurtzite that leads to flat bilayers because of 
a larger u parameter and an additional layer-parallel mirror 
plane. Furthermore the lattice constant c is only somewhat 
larger than the in-plane nearest neighbor distance, so that h- 
MgO is essentially five fold coordinated. 

To determine the equilibrium lattice parameters the total en- 
ergy was calculated for different cell volumes, while the cell- 
shape and internal parameters were allowed to relax. Using 
the Murnaghan equation of state (EOS) 11221 we obtained the 
energy- volume dependences E =E{V), the corresponding fits 
are represented in Fig.^ For the polymorphs of MgO, ZnO, 
and CdO being most stable in certain volume ranges around 
the equilibrium volumes, these fits lead to the equilibrium val- 
ues for the volume Vq and the total energy Eq per cation-anion 
pair, as well as the isothermal bulk modulus B, and its pres- 
sure coefficient B' = {dB/dp)p=o- The binding energy Eb has 
been calculated as the corresponding total energy Eq at zero 
temperature reduced by atomic total energies computed with 
spin polarization. In TableUthese parameters are summarized 
together with the lattice parameters we obtained. 

The E{V) curves in Fig.^clearly show that under ambient 
conditions the group-II oxides crystallize either in rs (MgO 
and CdO) or w (ZnO) structure. However, the binding en- 
ergies of the atoms in rs and w structure are rather similar, 
especially for CdO. The energy gains due to electrostatic at- 
traction on smaller distances (NaCl) and due to the better over- 
lap of sp^ hybrids (wurtzite) result in a sensitive energy bal- 
ance. Therefore we cannot give a simple explanation, why 
one of these crystal structures has to be favored over the other 
one. Furthermore we observe an anomalous structural trend 
along the cation row Mg, Zn, and Cd, which follows the 
anomalous trend of their covalent radii 1.36, 1.25, and 1.48 A 
]36]. Together with the oxygen radius of 0.73 A the result- 
ing nearest-neighbor distances 2.09, 1.98, and 2.21 A in the 
tetrahedrally coordinated wurtzite structure take a minimum 
for ZnO. That means besides the strong covalent bonds due 
to sp^ hybrid overlapping, also significant energy gain due 
to the Madelung energy occurs. As a result of the ab ini- 
tio calculations (cf. Table U the nearest-neighbor distances 
in the sixfold-coordinated rocksalt structure show a minor 
monotonous variation 2.13, 2.17, and 2.39 A. Also the cor- 
responding ionic energies, the repulsive interaction and the 
Madelung energy follow a monotonous trend. The sequence 
of the ratios of the nearest-neighbor distances of these two 
polymorphs with about 0.98, 0.91, and 0.93 may also be con- 
sidered as an indication why ZnO exhibits another equilibrium 
structure as MgO and CdO and, hence, for the nonexistence 
of an isostructural series. Both, the favorization of w or rs 
structure, as well as the cationic trend we observed, may be 
explained with the ideas of Zunger 13711 . Within his model 
he examined over 500 different compounds and predicted the 
correct equilibrium structures for the three materials we inves- 
tigated. 

Comparing the results with data of recent measurements or 
other first-principles calculations, we find excellent agreement 



TABLE I: Ground-state properties of equilibrium and high-pressure phases of MgO, ZnO, and CdO. All the experimental binding energies *) 
are taken from Ref. 1 23] as heat of vaporization or heat of atomization. 



Oxide 


Phase 




Lattice parameter 
(A or dimensionless) 


B (GPa) 


B' 


Binding energy 
(eV/pair) 


Reference 


MgO 


wurtzite 


a = 


c = u = 














3.322 


5.136 0.3916 


116.9 


2.7 


10.02 


this work 






3.169 


5.175 0.3750 


137 


- 


10.85 


theor. [201 




h-MgO 


a = 


c = u = 














3.523 


4.236 0.5002 


124.8 


4.3 


10.09 


this work 






3.426 


4.112 0.5000 


148 


- 


11.09 


theor. [22J 




rocksalt 


ao = 
















4.254 




148.6 


4.3 


10.17 


this work 






4.212 




- 


- 


10.26*) 


exp. [24J 






4.247 




169.1 


3.3 


10.05 


theor. [25J 






4.197 




169.0 


4.2 


- 


theor. [26J 






4.253 




150.6 


~ 


- 


theor. [27] 






4.145 




178.0 


- 


11.47 


theor. [22J 




cesium chloride 


ao = 
















2.661 




140.3 


4.1 


8.67 


this work 






2.656 




152.6 


3.4 


8.54 


theor. [^ 


ZnO 


wurtzite 


a = 


c = u = 














3.283 


5.309 0.3786 


131.5 


4.2 


7.20 


this work 






3.258 


5.220 0.382 


181 


4 


7.52*) 


exp. [22J 






3.250 


5.204 


183 


4 


- 


exp. [29] 






3.238 


5.232 0.380 


154 


4.3 


- 


theor. [2^ 






3.292 


5.292 0.3802 


133.7 


3.8 


7.69 


theor. [25J 






3.198 


5.167 0.379 


159.5 


4.5 


- 


theor. [30] 






3.183 


5.124 0.380 


162 


- 


10.64 


theor. [^ 




zinc blende 


flO = 
















4.627 




131.6 


3.3 


7.19 


this work 






4.633 




135.3 


3.7 


7.68 


theor. [2^ 






4.504 




160.8 


5.7 


- 


theor. [30] 




rocksalt 


ao = 
















4.334 




167.8 


5.3 


6.91 


this work 






4.275 




194 


4.8 


- 


exp. [32] 






4.287 




218 


4 


- 


exp. [28] 






4.271 




228 


4 


- 


exp. [29] 






4.283 




202.5 


3.5 


- 


exp. ^ 






4.272 




198 


4.6 


- 


theor. [28] 






4.345 




172.7 


3.7 


7.46 


theor. [2^ 






4.316 




175 


5.4 


- 


theor. [32J 






4.225 




209.1 


2.7 


- 


theor. [^ 






4.213 




210 


- 


10.43 


theor. [^ 




cesium chloride 


ao = 
















2.690 




162.4 


4.7 


5.76 


this work 






2.705 




156.9 


3.8 


6.33 


theor. [25J 


CdO 


wurtzite 


a = 


c = u = 














3.678 


5.825 0.3849 


92.7 


4.7 


5.97 


this work 






3.660 


5.856 0.3500 


86 


4.5 


5.30 


theor. [34] 




zinc blende 


ao = 
















5.148 




93.9 


5.0 


5.96 


this work 






5.150 




82 


3.0 


5.18 


theor. [34] 




rocksalt 


ao = 
















4.779 




130.5 


5.0 


6.00 


this work 






4.696 




148 


4 


6.40*) 


exp. [35] 






4.770 




130 


4.1 


5.30 


theor. [^ 
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FIG. 1: The nomalized total energy versus volume of one cation- 
oxygen pair. Several polymorplis liave been studied for MgO (a), 
ZnO (b), and CdO (c). w: thick sohd line, NiAs structure: thin solid 
line, /i-MgO structure: long dash-dotted line, rs: dashed line, zb: 
dotted line, CsCl structure: dash-dotted line 



(cf. Table|I|i. For ZnO the wurtzite ground-state and the NaCl 
(and CsCl) high-pressure phases are confirmed by experimen- 
tal studie s 123. [2^ BSl BSl B^ and other ah initio calcula- 
tions 1251123. l2y. 38, j91. Experimental results for the rocksalt 
ground-state of MgO |24] and CdO |27] exist, as well as other 
calculations for the equilibrium structure l25Ll26ll27l 13411 and 
the high-pressure phase (CsCl) |25, 34]. In the case of CdO 
also the CsCl structure has been studied experimentally 14011 . 
Our computed lattice parameters for the CdO polymorphs 
show excellent agreement, in particular compared with the re- 



sults of Jaffe et al. I25r . These authors also use a DFT-GGA 
scheme but expand the wave functions in localized orbitals of 
Gaussian form. However, also the agreement with values from 
other computations is excellent. Most of them use a DFT-LDA 
scheme which tends to an overbinding effect, i.e., too small 
lattice constants, too large bulk moduli, and too large binding 
energies. The overestimation of the binding energies of ZnO 
polymorphs in Ref. 1 3 1 ] is probably a consequence that not 
spin-polarized atomic energies have been substracted. 
The DFT-GGA scheme we used tends to underestimate 
slightly the bonding in the considered group-II oxide poly- 
morphs. This underestimation results in an overestimation of 
the lattice constants of about 1%. In the case of the c-lattice 
constants of w-ZnO and the ao constant of rs-QAO this dis- 
crepancy increases to roughly 2%. Also our calculated bulk 
moduli are always smaller than the experimental ones, which 
may be due to the mentioned underestimation. Besides the 
limitation of the computations we cannot exclude that sample- 
quality problems play a role for these discrepancies. For rs- 
MgO, w-ZnO, and r^-CdO the computed binding energies are 
close to the measured ones. The theoretical underestimation 
only amounts to 1, 5 or 7% which are small deviations. 



B. Pressure- and strain-induced pliase transitions 

Structural changes in the form of pressure-induced phase 
transitions are studied in detail in Fig. |2] for ZnO. Usually 
the Gibbs free energy G = f/ + pY — TS as the appropri- 
ate thermodynamic potential governs the crystal stability for 
given pressure and temperature. Its study however requires the 
knowledge of the full phonon spectrum. Therefore, we restrict 
ourselves to the discussion of the low-temperature limit, more 
strictly speaking to the electronic contribution to the enthalpy 
H ~ E + pV with the internal energy U{V)^E{V) and the ex- 
ternal pressure obtained as the trace of the stress tensor. The 
zero-point motional energy is neglected. Such an approach is 
sufficient for the discussion of the pressure-induced properties 
of relatively hard materials for temperatures below t hat g iven 
by the maximum frequency of the phonon spectrum 14111 . For 
a given pressure the crystallographic phase with the lowest 
enthalpy is the most stable one, and a crossing of two curves 
indicates a pressure-induced first-order phase transition. 
From Fig.|5]we derived the equilibrium transition pressure /?,. 
Using pt we obtained from p over V plots the initial vol- 
ume Vi and final volume Vf for the transitions, given here 
in units of the equilibrium volume Vo of the wurtzite poly- 
morph. We derive the values p, ^ ll.B. GPa (V,- = 0.92yo, 
Vf — 0.77Vb) for the transition wurtzite -^ NaCl and p, = 
261 GPa {Vi = O.SOVo, Vf = 0.47 Vo) for the transition NaCl -^ 
CsCl. The first values are in rough agreement with the experi- 
mental findings p, = 9.1 GPa (V, = 0.82Vo) ^33"], p, « 10 GPa 
128], or p, w 9 GPa |42], though our calculations indicate 
a slightly higher stability of the wurtzite structure over the 
rocksalt one. The computed pt value is in reasonable agree- 
ment with other calculations (see I2al and references therein). 
Another pressure-induced phase transition between NaCl and 
CsCl structure is found at a transition pressure of 261 GPa, 
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TABLE II: Lattice parameters of the energetically preferred poly- 
morphs of MgO and CdO under biaxial strain along the bilayer- 
stacking axis. 
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25 50 75 " 200 225 250 275 ^"^ 
Pressure p (GPa) 

FIG. 2: Enthalpy per cation-anion pair of ZnO phases as a function 
of hydrostatic pressure, w: thick solid line, rs: dashed line, CsCl 
structure: dash-dotted line. The curve for zinc blende is practically 
identical with that of wurtzite. 



very similar to a previous calculation 12511 which predicted a 
value of pt — 256 GPa. 

Applying the common-tangent method the E{V) curves in 
Figs.^ and^ for MgO and CdO already show that pressure- 
induced phase transitions are hardly observable or should oc- 
cur at high transition pressures (and hence small transition 
volumes). In the case of MgO the crossing of the enthalpies 
gives a value of p, ~ 508 GPa for the transition NaCl — > CsCl 
structure. For CdO we obtain pt w 85 GPa for the transition 
NaCl -^ CsCl structure. In comparison with the result of a 
previous calculation |27] we find good a gree ment. Another 
predicted value amounts to 515 GPa (see ll25ll and references 
therein). Also the CdO value is in excellent agreement with 
other DFT-GGA calculations p, ^ 89 GPaJpl and even with 
the result of measurements p, — 90 GPa |40]. 
For MgO (cf. Fig.^) we have to mention an interesting re- 
sult for negative pressures: at large volumes of about 1.3Vo, 
indeed wurtzite is more stable than the NaCl structure. How- 
ever, when we start the atomic relaxation with the wurtzite 
structure and slightly decrease the cell volume, we observe a 
transition into the /i-MgO structure. The corresponding total 
energy minimum therefore lies between the rocksalt and the 
wurtzite minima at a volume of about V = I .IVq. Around this 
volume there is no energy barrier between the wurtzite and 
h-MgO structures and the wurtzite geometry only represents 
a saddle point on the total energy surface, whereas we ob- 
serve h-MgO to be an intermediate structure on the way from 
wurtzite to rocksalt, as discussed in 12011 . There is another 
indication for this transition: A decrease of c/a leads to an 
increase ofw, followed by a sudden relaxation into the h-MgO 
structure 14311 . Comparing the u parameters of our wurtzite 
structures (cf. Table|IJ, we observe that w-MgO has the high- 
est u and so this relaxation is most probable for MgO. 

An important point for the above-mentioned heterostruc- 
tures it is the possibility to grow pseudomorphically one ma- 
terial on the other. Because wurtzite ZnO substrates are 
commercially available, the question arises if such growth of 
certain polymorphs of MgO or CdO on a w-ZnO substrate 



w-MgO 


r^-CdO 


a (A) 


c(A) 


u 


a (A) 


f(A) 


3.283 


5.201 


0.3856 


4.643 


8.512 



with [0001] orientation is possible. For that reason we com- 
pare the fl-lattice constant of w-ZnO, a = 3.283 A, with the 
corresponding lattice constants a of hexagonal modifications 
of MgO or CdO and the second-nearest neighbor distances 
flo/V2 in the cubic cases. The corresponding values are 3.523 
(/z-MgO), 3.322 (w), and 3.008 A (rs) for MgO or 3.678 (w), 
3.640 izb), and 3.379 A {rs) for CdO. With these values the 
resulting lattice misfits are 7.3, 1.2, and —8.4% for MgO and 
12.0, 10.9, and 2.9% for CdO for [0001]/[111] interfaces. 
From the point of lattice-constant matching pseudomorphic 
growth of w-MgO and, perhaps, also r^-CdO should be pos- 
sible on a w-ZnO substrate. This conclusion is confirmed by 
the total-energy studies (see Fig.|3l for several MgO and CdO 
polymorphs biaxially strained in [0001] or [111] direction of 
cation-anion bilayer stacking. For the curves in Fig.|3lwe kept 
the corresponding a-lattice constant fixed at the w-ZnO value 
and computed the total energy per cation-anion pair for several 
values of the c-lattice constant of the coiTesponding hexagonal 
crystal or of the resulting rhombohedral crystal, the resulting 
lattice parameters are given in Tablellll 

The results for MgO are most interesting: In the presence 
of a weak biaxial strain in [0001] direction of about 1.3% the 
most energetically favorable geometry is the wurtzite struc- 
ture with a resulting c-lattice constant of c = 5.21 A. The two 
polymporphs also considered here, the /z-MgO structure and 
the trigonally distorted rs geometry, are much higher in en- 
ergy. Furthermore the energetical ordering of the MgO poly- 
morphs (cf. Fig.^) is completely changed and MgO adopts 
the wurtzite crystal structure of the substrate which is differ- 
ent from its rocksalt geometry in thermal equilibrium. From 
the u parameters for constrained w-MgO (see Table |llji and 
w-MgO in equilibrium (see Table |I} we find that according 
to 1 43] the biaxially strained wurtzite structure should be more 
stable against a transition into /i-MgO, because its u parameter 
is closer to its ideal value. We conclude that pseudomorphic 
growth of MgO should be possible on ZnO(OOOl) substrates. 
We find a somewhat different situation for CdO. Apart from 
the large lattice misfit which probably prohibits pseudomor- 
phic growth, its lowest energy structure is still derived from 
the rs atomic aiTangement. Of course, there is a strong trigo- 
nal distortion giving rise to a c-lattice constant of c = 8.512 A 
in comparison with the thickness \/3flo of three CdO bilayers 
in [111] direction at thermal equilibrium with c = 8.277 A. 
In the case of CdO the two other crystallographic structures, 
w and zb, are energetically less favorable also in the biaxially 
strained case. 
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FIG. 3: Total energy of biaxially strained MgO (a) and CdO (b) 
polymorphs versus the c-lattice constant in [111] or [0001] direction 
for an a-lattice constant fixed at the value a = 3.283 A of wurtzite 
ZnO. w. thick solid line, rs: dashed line, zb: dotted line, and h-MgO 
structure: long dash-dotted line. 



IV. EXCITED-STATE PROPERTIES 
A. Band structures and electronic densities of states 

To study the influence of the atomic geometry, more pre- 
cisely the polymorph of the group-II oxide compound, we 
show in Figs. |4] |5j and |6] the band structures as calculated 
within DFT-GGA for the most stable polymorphs in a certain 
volume range around the equilibrium volume (cf. Sec. III). 
These are the rocksalt and cesium chloride structures, supple- 
mented by the wurtzite structure and the zinc -blende structure. 
By examining the electronic structure in the KS approach one 
neglects the excitation aspect 11441 mil and, hence, underesti- 
mates the resulting energy gaps and interband transition en- 
ergies. In GW approximation the corresponding quasiparticle 
(QP) ener gy co rrections due to the exchange-correlation self- 
energy f44','45','4^ amount to about 2.0 eV for ZnO [47] and 
3.6 eV for MgO |48]. Nevertheless the independent-particle 
approximation |49] is frequently a good starting point for the 
description of optical properties. 

For MgO the band structures and densities of states in Fig.|3 
indicate an insulator or wide-gap semiconductor independent 
of the polymorph. Except the CsCl-structure, which has an 
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FIG. 4: Band structure and density of states (normalized per pair) 
for MgO polymorphs calculated within the DFT-GGA framework: 
(a) wurtzite, (b) zinc blende, (c) rocksalt, and (d) cesium chloride. 
The shaded region indicates the fundamental gap. The valence band 
maximum is chosen as energy zero. 



indirect gap between M and F, they all have direct fundamen- 
tal gaps at the F point in the BZ. We find rather similar values 
4.5 eV (rs), 3.5 eV (zb), 4.2 eV (NiAs), and 3.3 eV (/z-MgO) 
for the gaps of cubic and hexagonal crystals despite the dif- 
ferent lattice constants, coordination, and bonding. There are 
also similarities in the atomic origin of the bands. The 02s 
states give rise to weakly dispersive bands 15-18 eV below 
the valence-band maximum (VBM). They are separated by 
an ionic gap of 10-12 eV from the uppermost valence bands 
with band widths below 5 eV. Due to the high ionicity of the 
bonds the corresponding eigenstates predominately possess 
02p character. For the same reason the lowest conduction 
bands can be traced back to Mg3s states. Also mentionable is 
the relatively weak dispersion of the uppermost valence bands 
for the wurtzite and zinc -blende structures which leads, com- 
pared with the rocksalt and CsCl structures, to a quite high 
DOS near the VBM. 

For the equilibrium rocksalt polymorph we calculate besides 
the direct gap at F other direct gaps at X and L in the BZ: 
Eg{X) = 10.43 eV and Ef,{L) = 8.37 eV. These values are 
in good agreement with results of other DFT calculations us- 
ing an LDA exchange-correlation functional and a smaller lat- 
tice constant |48, 50, 51]. In their works |48|/|50] these au- 
thors computed QP gap openings of about 3.59/2.5 eV (F), 
3.96/2.5 eV (X), and 4.00/2.5 eV (L) in rough agreement 
with the values 3.06 eV of Shirley [51] and 2.94 eV derived 
from the simple Bechstedt-DelSole formula for tetrahedrally 
bonded crystals ]52]. With these values one obtains QP gaps 
which approach the experimental values of ^.^(F) —1.1 eV, 
£UX) = 13.3 eV, and £g(L) = 10.8 eV 16] or4(r) =7.83eV 
II53I1 . A more sophisticated QP value for the fundamental gap 
is£g(F)=7.79eV154]. 

The electronic structures of ZnO plotted in Fig.|5]show sev- 
eral similarities to those of MgO. However, the 02s bands 
now appear at about 15.5-18.5 eV below the VBM, and the 
uppermost valence bands of predominantly 02p character are 
found in the range from to —4 eV The lowest conduction- 
band states (at least near F) are dominated by Zn4s states. In 
comparison with the fundamental energy gaps of MgO those 
of ZnO are smaller We compute 0.73 eV (w), 0.64 eV (zb), 
1.97 eV (rs) for the direct gap at F. However, for the rocksalt 
geometry the VBM occurs at the L point and therefore this 
high-pressure phase is an indirect semiconductor with a gap 
of Ef, (F — L) = 0.75 e V. This observation is in agreement with 
findings by room-temperature absorption measurements and 
DFT-LDA calculations |25l |47l |5j. Another local valence 
band maximum which is almost as high as the one at L occurs 
at the E line between K and F. 

New features, that are not observable for MgO are caused by 
the Zn3d states. These shallow core states give rise to two 
groups of bands (at F) clearly visible in the energy range of 
4-6 eV below the VBM, which generally show a splitting and 
a wave-vector dispersion outside F. The huge peaks caused 
by these basically Zn3£/-derived bands are clearly visible in 
the DOS. Furthermore, the Zn3ii states act more subtle on 
the band structure via the repulsion of p and d bands caused 
by the hybridization of the respective states. The effects of 
this pd repulsion can be discussed most easily for the F-point 
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FIG. 5: Band structure and density of states (normalized per pair) 
for ZnO polymorphs calculated within the DFT-GGA framework: 
(a) wurtzite, (b) zinc blende, (c) rocksalt, and (d) cesium chloride. 
The shaded region indicates the fundamental gap. The valence band 
maximum is chosen as energy zero. 



in the BZ: In the case of the zb polymorph (with tetrahedral 
coordination and tJ symmetry), the hybridized anion p{t2) 
and cation d{t2) levels give rise to the threefold degenerate 
ri5{pd) and ri5{dp) levels po). In the case of wurtzite the 
levels are doubled at F (with respect to the zfe-polymorph) due 
to the band folding along the [111]/[0001] direction. In ad- 
dition, these states are influenced by a crystal-field splitting. 
Moreover, in the case of the tetrahedrally coordinated zb and 
w polymorphs the pd repulsion reduces the fundamental di- 
rect gaps at r as the ri5{pd) are pushed to higher energies, 
while the Fi^ conduction band minimum (CBM) remains un- 
affected. In the case of the rs polymorph (with sixfold coor- 
dination and octahedral O^^ symmetry), the anion p{ti) level 
gives rise to the threefold degenerate Fi5(p) valence band, 
while the symmetry-adapted e (twofold degenerate) and t2 
(threefold degenerate) combination of the cation d states gen- 
erates the Fi2(c/) and F25'(ii) bands. The respective states do 
not hybridize and the pd repulsion as well as the correspond- 
ing gap shrinkage vanishes |56]. However in other regions of 
the BZ, the bands are subject to the pd repulsion and thereby 
raised at points away from F. Consequently the pd repulsion, 
or more exactly its symmetry forbiddance at F is the reason 
why the rs polymorph is an indirect semiconductor. Com- 
pared to InN the pd repulsion is larger but does not give rise 
to a negative Fi {s)-ri5{pd) gap as found there 15711 . 
Using QP corrections the small energy gaps obtained within 
DFT-GGA are significantly opened. With appropriate param- 
eters as bond polarizability 0.78, nearest-neighbor distance 
2.01 A, and electronic dielectric constant 4.0 |58]) for the 
tetrahedrally coordinated ZnO the Bechstedt-DelSole formula 
gives a QP shift of about 1 .95 eV. This gap correction is brack- 
eted by other approximate values of 0.85 eV |59], 1.04 eV 
ll60ll . and 2.49 eV |61]. However, it agrees well with the 
value 1.67 eV derived within a sophisticated QP calculation 
ll47ll . Our values resulting for the direct fundamental gap at F 
2.68 eV (w) or 2.59 eV (zb) are larger than the QP gap value of 
2.44 eV f^Tl but still clearly underestimate the experimental 
gap of 3.44 eV [62, 63]. This underestimation is sometimes 
related to the overscreening within the random-phase approx- 
imation (RPA) used in the QP approach 147). We claim that 
one important reason is the overestimation of the pd repulsion 
and, hence, the too high position of the VBM in energy due to 
the too shallow d bands in LDA/GGA. Even when including 
many-body effects, the Zn3d bands are still to high in energy. 
For a better treatment of the QP effects also for the semicore d 
states, the pd repulsion should be reduced, which requires the 
inclusion of non-diagonal elements of the self-energy operator 
or the use of a starting point different from LDA/GGA, e.g. a 
generalized KS scheme [64]. A recent work [55[ reported that 
rocksalt ZnO has an indirect band gap of 2.45±0.15 eV mea- 
sured from optical absorption. Using the above estimated QP 
shift of 1.95 eV and the indirect gap in DFT-LDA quality of 
0.75 eV one finds a QP value of 2.60 eV close to the measured 
absorption edge. 

In the case of the same polymorph the band structure and 
density of states of CdO in Fig. |6] shows several similarities 
with those for ZnO. At about 15.5-17 eV below the VBM 
occur the 02s bands, and the Cd4ii bands are observed in 
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FIG. 6: Band structure and density of states (normalized per pair) 
for CdO polymorphs calculated within the DFT-GGA framework: 
(a) wurtzite, (b) zinc blende, (c) rocksalt, and (d) cesium chloride. 
The shaded region indicates the fundamental gap. The valence band 
maximum is chosen as energy zero. 



the energy interval from —5.9 to —6.6 eV. The uppermost 
02p-derived valence bands possess a maximum band width 
of about 3.3 eV, a value smaller than that found in recent pho- 
toemission studies |65]. Usually the conduction bands are 
well separated from the valence bands. However, for all poly- 
morphs the lowest conduction band shows a pronounced min- 
imum at the center of the BZ. Within an energy range of a 
few tenths of an electron volt above its bottom, this band is 
isotropic but highly nonparabolic. 

In rs-CdO we find a direct gap of about 0.66 eV at the F 
point. As in the case of rs-ZnO also for r^-CdO the max- 
ima of the valence bands occur at the L point and at the E line 
between F and K. These maxima lie above the CBM, which 
results in negative indirect gaps of about —0.5 1 and —0.43 eV, 
and therefore our band structure indicates a halfmetal. As de- 
scribed for rs-ZnO, the pd repulsion is responsible for this ef- 
fect ll66ll . Our band structure is in qualitative agreement with 
a DFT-LDA calculation f6y], in particular with respect to the 
band dispersions. However, Ref. |65] found positive direct 
and indirect gap values. The gap opening of about 1 eV with 
respect to our values may be a consequence of the used basis 
set restricted to a few Gaussians. Corresponding experimental 
values are 0.84 and 1.09 eV for the indirect gaps and 2.28 eV 
for the lowest direct gap at F. Their comparison with the val- 
ues calculated within DFT-GGA indicates effective QP gap 
openings of about 1.3-1.7 eV. However, these values should 
be considerably influenced by the pd repulsion, at least out- 
side the BZ center. While the Cd4d bands are about 6.4 or 
6.6 eV below the valence band maximum at F, the experi- 
mental distances are about 12.4 or 13.3 eV (with respect to 
the Fermi level) ll67ll . A more recent measurement indicates 
an average binding energy of Cd4d relative to VBM of about 
9.4 eV |65]. Altogether the band structures of CdO are rather 
similar to those of InN |68|. Within DFT-GGA in both cases 
small negative s-p fundamental gaps are found near F. The 
main difference is related with the position of the 4d bands, 
because the binding energy of In4ii electrons is larger than 
that of Cd4d electrons. 

The band structures presented in Figs. |4] |5j and |6l show 
clear chemical trends along the series MgO, ZnO, and CdO 
for a fixed crystal structure. To clarify these trends we list 
in Table ^^ some energy positions for the two most impor- 
tant polymorphs - rocksalt and wurtzite. Among these ener- 
gies we also list the lowest conduction and highest valence 
bands for two high-symmetry points in the corresponding fee 
(F, L) or hexagonal (F, A) BZ. In the case of the cubic sys- 
tems the position of the twofold (threefold) degenerate shal- 
low core d level F12 (F25/) at F is given. For the rocksalt 
phase as well as the wurtzite polymorph the level positions 
follow a clear chemical trend with the cations Mg, Zn, and 
Cd. The energetical position of the empty conduction-band 
levels as well as the filled valence-band states decreases with 
respect to the VBM. Consequently, the average gaps decrease 
along the series MgO, ZnO, and CdO. At least for tetrahe- 
drally coordinated compounds the general trend is governed 
by both, the splitting of the s- and p-valence energies as well 
as the nearest-neighbor distances |58]. In the case of wurtzite 
the effect of the non-monotonous behaviour of both quantities 



TABLE III: Characteristic energy levels (in eV) in the band structure 
of the rocksalt and wurtzite polymorphs of MgO, ZnO, and CdO as 
calculated within the DFT-GGA framework. The center F of the BZ 
and an L and A point at the BZ surface in [111]/[0001] direction 
are chosen for the Bloch wave vectors. The lowest conduction (f) 
bands and the highest valence (v) bands are studied. In the case of 
rocksalt also the positions of the Fn and r25' d bands are given. The 
uppermost Fis,. or T(,y valence band at F is used as energy zero. 
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may give a monotonous net effect. Another atomic tendency, 
the increasing cation-/?-anion-(i splitting from ZnO to CdO 
I5a.l5«l . is correlated with the increase of the distance of the 
d bands to the VBM. 

We have also computed the volume deformation potentials for 
the direct gaps at F in rocksalt MgO to ay = —9.39 eV and in 
wurtzite ZnO to ay = — 1 .55 e V, as well as for the indirect gap 
between L and F in rocksalt CdO to ay = —1 .87 eV The abso- 
lute value of the gap deformation potential of MgO is clearly 
larger than those of ZnO or CdO, due to the stronger bonding 
in MgO. Despite different gap states and ground-state poly- 
morphs the values for ZnO and CdO are rather similar. Their 
smallness is a consequence of the large ionicity and the rela- 
tively large bond length. According to the estimate of 1 8], this 
may explain the small gaps of ZnO and CdO in comparison 
to ZnS and CdS. The deformation potential ay = —1.55 eV 
calculated for w-ZnO within DFT-GGA is much smaller com- 
pared to the experimental value of ay = —3.51 eV |69]. The 
reason should be the neglect of the large QP corrections not 
taken into account. 



B. Dielectric functions and optical properties 

For computing the dielectric function £(©) we have cho- 
sen the wurtzite and rocksalt polymorphs since they give the 
equilibrium geometries. In addition, we show spectra for 
zinc -blende crystals because of the mentioned similarity with 
wurtzite, but without its optical anisotropy. In Figs.0|8j and|9] 
the influence of the crystallographic structure and of the anion 
on the dielectric function of a group-II oxide is demonstrated. 
We observe that the change of the coordination of the atoms 
has a strong influence, which perhaps is best noticeable in 
the imaginary parts of the dielectric function. Going from 
the fourfold (w, zb) to the sixfold (rs) coordination, the ab- 
sorption edge shifts towards higher photon energies and the 
oscillator strength increases. Thereby the screening sum rule 
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is less influenced. One observes a clear increase of the (high- 
frequency) electronic dielectric constant £„ = 5R (e(0)) along 
the series MgO, ZnO, and CdO, so the main influence is due 
to the chemistry. For the rs structure we compute £00 = 3.16, 
5.32, 7.20 in qualitative agreement with the reduction of the 
fundamental gap. Thereby, for CdO the accuracy is reduced 
due to the difficulties with the DFT-GGA band structure dis- 
cussed above. 

In the real parts there are several features that are similar, 
independent of the polymorph. There is only small variation 
of the £00 with the different polymorphs. In the case of ZnO 
we find £oox,v = 5.24, £«,jj = 5.26 (w), £„ = 5.54 (zb), and £„ = 
5.32 (rs). These values are larger than the constants £ooxv = 
3.70 and £00,,; = 3.75 measured for w-ZnO |7]. This fact may 
be traced back to the underestimation of the fundamental gap 
within DFT-GGA. With the values £„ = 2.94 (measured iJl) 
and £00 = 3.16 (this work) the agreement is much better for 
rs-MgO. 

The lineshapes resulting from the dielectric function are 
discussed in detail for ZnO (Fig.|8|i. For smaller frequencies 
the curves 3i{e{co)) exhibit maxima close to the absorption 
edge. These maxima are followed by regions with the general 
tendency for reduced intensity but modulated by peak struc- 
tures related to critical points in the BZ. For frequencies larger 
than about 12.5 eV the real part becomes negative for all poly- 
morphs. For the imaginary part in Fig.|8]the differences are 
larger - in agreement with the results of other calculations 
II39I1 . In the case of rocksalt there is a monotonous increase to 
the first main peaks at hco — 5.5 and 7.0 eV which should be 
shifted towards higher energies in the experimental spectra ac- 
cording to the huge QP shifts discussed for the band structure. 
These two peaks are due to transitions near the L and X points 
of the fee BZ (cf. Fig.|5}. Therefore they can be classified as 
El and E2 transitions |70]. The subsequent peaks at about 9.5, 
12.0, and 13.2eV should be related to iij (i.e., second valence- 
band into lowest conduction-band at L) and £2 transitions. 
We find a different situation for the wurtzite and zinc-blende 
structures. One observes a steep onset of the absorption just 
for photon energies only slightly larger than the fundamental 
band gap. In a range of about 4 eV it follows a more or less 
constant or even concave region. Such a lineshape is clearly 
a consequence of the pronounced conduction-band minimum 
near F, the nonparabolicity of the conduction band and the 
light-hole valence band. It has also been observed experimen- 
tally iTlll . We therefore have a rather similar situation in the 
case of InN |68!l, the III-V compound with constituents neigh- 
boring CdO in the periodic table of elements. 
Within a four-band k • p Kane model one finds for the imagi- 
nary part of the dielectric function in the case of cubic systems 

m 
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FIG. 7: Real and imaginary part of the frequency-dependent dielec- 
tric function for MgO polymorphs wurtzite (a), zinc blende (b), and 
rocksalt (c) as calculated within the independent-particle approxima- 
tion using Kohn-Sham eigenstates and eigenvalues from DFT-GGA. 
Imaginary part: solid line, real part: dotted line. In the case of 
wurtzite besides the tensor components Evx(ft)) = £yy{(0) also the 
zz-component Ezzico) is presented (imaginary part: dashed line, real 
part: dash-dotted line). 
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with the Bohr radius qb, the fundamental direct band gap 
Eg, and the characteristic energy Ep, which is related to the 
square of the momentum-operator matrix element between s 
and p valence states. Indeed formula Q leads to a constant 



3 (e(«)) = 3(e^/2flB£'p)2 for hco > Eg, i.e., away from the 
absorption edge. Interestingly the absolute plateau values in 
Figs.|8t andlsj) of about 1.8 may be related to an energy E^ 
of about w 37 eV, a value which is not too far from those in 
other estimations |72, 73]. The same k • p model gives for the 
conduction-band mass near F the expression m* = mo/[l + 
Ep/Eg]. With an experimental gap energy of Eg ^ 3.4- eV 
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FIG. 8: Real and imaginary part of the frequency-dependent dielec- 
tric function for ZnO polymorphs wurtzite (a), zinc blende (b), and 
rocksalt (c) as calculated within the independent-particle approxima- 
tion using Kohn-Sham eigenstates and eigenvalues from DFT-GGA. 
Imaginary part: solid line, real part: dotted line. In the case of 
wurtzite besides the tensor components exi-(tt)) = eyv(ft)) also the 
zz-component Z^jioi) is presented (imaginary part: dashed line, real 
part: dash-dotted line). 



FIG. 9: Real and imaginary part of the frequency-dependent dielec- 
tric function for CdO polymorphs wurtzite (a), zinc blende (b), and 
rocksalt (c) as calculated within the independent-particle approxima- 
tion using Kohn-Sham eigenstates and eigenvalues from DFT-GGA. 
Imaginary part: solid line, real part: dotted line. In the case of 
wurtzite besides the tensor components Exi^®) = evv(^) also the 
zz-component ZTjjyOi) is presented (imaginary part: dashed line, real 
part: dash-dotted line). 



it results an electron mass of about m* = 0.09mo somewhat 
smaller than the experimental value of about 0.1 9mo |73] or 
0.28mo |7]. From the conduction band minima plotted in 
Figs. |5^ and|5j3 we derive band masses in DFT-GGA qual- 
ity of about m* = 0.153mo (zb) and m* = O.lSlmo/O.lSOmo 
(w) with a negligible anisotropy. Absorption peaks occur at 
higher photon energies at about 7.5, 9.1, and 10.6 eV for zb- 



ZnO and 6.5, 9.9, and 10.7 eV for w-ZnO. These peaks may 
be related to E1/E2, E[, and £2 transitions in the case of the 
zinc-blende structure. For the wurtzite polymorph the inter- 
pretation is more difficult, however, the lowest peak should 
mainly be caused by transitions on the LM line between the 
highest valence band and the lowest conduction band. 

We have to mention that the influence of many -body effects 
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such as QP shifts |44, 45, 46] and exci tonic effects I74ll has 
been neglected. The QP effects lead to a noticeable blueshift 
of the absorption spectra, while the exci tonic effects, basically 
the electron-hole pair attraction, give rise to a redshift and a 
more or less strong mixing of interband transitions. Further- 
more, the neglected effects should cause a redistribution of 
spectral strength from higher to lower photon energies I68il . 



V. SUMMARY AND CONCLUSIONS 

Using the ab initio density functional theory together with 
a generalized gradient corrected exchange-correlation func- 
tional, we have calculated the ground-state and excited-state 
properties of several polymorphs of the group-11 oxides MgO, 
ZnO, and CdO. We have especially studied the rocksalt and 
wurtzite structures which give rise to the equilibrium geome- 
tries of the oxides. Zinc blende has been investigated to study 
the same local tetrahedron bonding geometry but without the 
resulting macroscopic anisotropy. Cesium chloride is an im- 
portant high-pressure polymorph. In addition, two hexagonal 
structures, nickel arsenide and /i-MgO, have been studied for 
MgO. 

In agreement with experimental and other theoretical find- 
ings the rocksalt (wurtzite) structure has been identified as 
the equilibrium geometry of MgO and CdO (ZnO). The non- 
isostructural series of the three compounds with a common 
anion has been related to the non-monotonous variation of 
the cation size along the column Mg, Zn, and Cd. We cal- 
culated binding energies which are in good agreement with 
measured values. MgO and CdO undergo a pressure-induced 
phase transition from the NaCl into the CsCl structure. w-ZnO 
first transforms into rocksalt geometry but also shows a phase 
transition into the cesium chloride geometry at higher hydro- 
static pressures. Our values for the transition pressures agree 
well with predictions of other ab initio calculations and exper- 



imental observations. We predicted the possibility of pseudo- 
morphic growth of MgO in biaxially strained wurtzite geom- 
etry on a w-ZnO substrate. 

The atomic coordination and hence the polymorph has a 
strong influence on the distribution of the allowed Bloch ener- 
gies. This fact has been clearly demonstrated by the compar- 
ison of band structures and densities of states computed for 
different crystallographic structures. This fact holds in partic- 
ular for the fundamental energy gaps which however strongly 
suffer from the neglect of the excitation aspect, the so-called 
QP corrections. In the case of ZnO and CdO the semicore 
Zn3d and Cd4d states also contribute to the gap shrinkage. In 
the framework of DFT-GGA they are too shallow and hence 
give rise to an overestimation of the pd repulsion (which shifts 
the uppermost p-like valence band towards higher energies). 
The effect of the crystal structure on the dielectric function is 
weaker as in the case of the band structures due to the occur- 
ring Brillouin-zone integration. We found a stronger influence 
of the cation. A clear chemical trend has been observed for the 
electronic dielectric constants. Along the series MgO, ZnO, 
and CdO also the averaged spectral strength increases. The 
crystal structure has the most influence on the lineshape of the 
absorption edge. This effect has been intensively discussed 
for ZnO. For the wurtzite and zinc -blende polymorphs we ob- 
serve a steep onset in the absorption, followed by a plateau- 
like frequency region. 
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